arXiv: 1505.06509vl [math.CA] 16 May 2015 


Exact analytical approach to differential equations 
with variable coefficients 


Mauro Bologna 

Institute de Alta Investigacion, Universidad de Tarapaca-Casilla 6-D Arica, Chile 
E-mail: mauroh69@gmail.com 

Abstract. This paper shows how to build a formal analytical solution for a 
differential equation of arbitrary order and with variable coefficients. It proofs that 
the most known approximated solutions for such a problem can be derived from the 
analytical expression presented in the paper. The formalism can be easily extended 
to the infinite dimensional case such as the quantum time-dependent Hamiltonian 
problem. 


TAGS numbers: 02.30.Hq, 03.65.Fd 


Exact analytical approach to differential equations with variable coefficients 

1 . Introduction 


2 


Linear differential equations are ubiquitous in all areas of science. Very typically, a 
scientific problem is described by a mathematical formalism and, very likely, eventually 
this implies a differential equation, possibly linear, that is to say 

+ an-i{z)y^"‘~^\z) H-h ao{z)y{z) = 0 (1) 

where y^^'^ means the jth derivative of y{z). It is virtually impossible to list a 
complete bibliography on the topic. Its applicability ranges from quantum mechanics [1], 
astrophysics [2] , chemistry and stochastic processes [31111 El E] etc. just to cite text books 
or paper collection where the reader can find problems ending into linear differential 
equation with variable coefficients. Also integro-differential and non-linear differential 
equations can be reduced to a finite or inhnite series of linear differential equations. 
A very famous example is the Riccati equation, s'{z) + s^z)"^ = f{z), that, via the 
transformation = y'{z)/y{z), can be reduced to the second order differential 

equation y"{z) = f{z)y{z). 

It is well known that Eq. ([T]) can be written as a hrst order equation in the matrix 
formalism, i.e. 

f-Yiz) = Y(z)M(z) ( 2 ) 

dz 

or, alternatively, 

fY(z) = M{z)Y(z) ( 3 ) 

dz 

where Y(z) is the vector Y(z) = {y{z), • • •, y^^~^\z)}. More in general the matrix M{z) 
describes a linear system of differential equations and may have an inhnite dimension 
such as in the time dependent Hamiltonian problem. Indeed Eqs. ([2]) and (j3]) play a 
fundamental role in physics when the matrix M(z) is the Hamiltonian of the system. In 
particular Eqs. (Ej) and Eq. (|3]) may describe either the stationary Schrodinger equation 
or the time-dependent Schrodinger equation, i.e. the two equations 

( 9 ^ 

- + V{x)'ijj{x) = E'ijj{x) (4) 

2m ox^ 

t) = H{t)f){x, t) (5) 

The aim of this paper is to hnd a formal analytical solution of Eq. ([2]) [or Eq. (Ej)]. It will 
be shown that, starting from the analytical solution, it is possible to recover and improve 
the most important approximations known in literature as part of one mathematical 
scheme (see Ref. mi as excellent review and text books about approximated method). 

The paper is organized as follows. Section [2] shows how to build a formal solution 
for a differential equation with variable coefficients. Section El presents few example of 
exact solutions obtained using the methods of Section El Section E] applies the results of 
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the previous sections to the time dependent Hamiltonian problem. Section E] shows how 
to recover the solution of an ordinary differential equation using the time-dependent 
Hamiltonian approach previously developed. Section |6] shows and tests an analytical 
formula (derived from the previous sections) to evaluate the eigenvalues in the problem 
of a differential equation with boundary conditions. Finally Section [7] summarizes the 
obtained results and the future work. 


2. Formal solution for linear n-order differential equations with variable 
coefficients 


To solve the problem set by Eq. da. we may solve the equivalent problem under the 
form of Eq. ([2|) [or Eq. ([3])]. Our starting point is to evaluate the high order derivatives 
of Eq. ([2]). Considering the second derivative we have 

Y^^\z) = Y{z)M{z) + Y{z)M‘^{z) = Y{z)[OM{z)] (6) 


where we defined 

0 = ^ + M{z). 
dz 

It is straightforward to show that, for n > 1 


yO+i)(^) _ y{z)[0'^M{z)]. 


(7) 

( 8 ) 


Performing the Taylor expansion of the vectorial function Y[zo + z) at the point z = zo, 
we may write 


Y{zo + z) = J2y^^\^o)^ 

n=0 


Y(zo) + V(zo) 

n=0 


Qn 

{n + 1)! 


M{zo) 




Without loss of generality we set = 0 and, after few algebraic steps, we can write the 
formal solution for Y [z) as 


Y(z) = Y(Q) 


1+ (ixp\uO]M{zo) du 


( 10 ) 


L ^0 J 

where I is the identity matrix and it is understood that O acts on the right, i.e. on M{zq). 
The symbol exp [uO] M[zf) |^g=o means that it must be evaluated at zq = 0. Analogously 
we may write the solution for 


fL 

dz 


Y(z) = M{z)Y{z) 


(11) 


as 


Yiz) 


1 + 


M{zo)exp [uO] 1^0=0 du 


r(o) 


( 12 ) 


L ^'O J 

where it is understood that O acts on the left. Being the operator O the sum of two 
operators we can use several tools present in literature to manipulate the solution given 
by Eqs. (nnii or (na. Mostly there are two important ways to proceed: the Baker- 
Campbell-Hausdorff formula [HI [101 [H] (and its dual, the Zassenhaus formula [12]) and 
the Trotter product formula [13] and consequently the path integral technique [H] . The 
former is left to a future work while, in this paper, we will focus on the latter. 
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3. Linear n-order differential eqnations with variable coefficients: exact 
examples 


In this section we consider examples where the formal solution gives exact results. 
As first example let us consider the case when M{z) = Mq i.e. a constant matrix. 
Since [^,Mo] = 0, we may factorize the exponential 



' f d 


d' 

exp 


Mq = exp 



exp [uMf\ Mo = exp [uMf\ Mq, 
and resolving the integral of Eq. (ITOll . we obtain the well known result 


Y{z) = F(0) exp [zMq] . 


(13) 


Let us focus on the case where M{z) has polynomial elements. Polynomial coefficients in 
differential equations are widely studied in literature and there are many works on this 
topic (here a non exhaustive list of references [151 [HI [HI [HI [IS [2S [211 [22] )• We limit 
ourselves to few exact example in the case when the matrix elements are polynomials 
satisfying to certain requirements. Setting constraints on the matrix M{z) we can find 
analytical expression for Eq. (|T0|1 . For example, requiring that 


M{zf 


constant/, 


dz'^ 


M{z) = 0, 


and using the properties of the operator O, it can be showed that 


(14) 


O^^M{z) = M{z){OM{z)Y, (15) 

0^^+^M{z) = {OM{z)Y+^. (16) 


Defining the matrix A‘^{z) = OM{z) we obtain 

du + Aq / sinh 

Jo 

where for compactness Mq = M(0) and Aq = A(0). If Aq is an invertible matrix, then 
expression (HZD may be written 

Y{z) = Y (0) [MqAq ^ sinh [zAq] + cosh [zAq]\ . (18) 

Note that we did not specify the dimension of the matrix M{z). This implies that 
solution flTTl) holds for an arbitrary order of differential equations described by conditions 
O- More: the second condition fflTl) implies that the elements of the matrix are hrst- 
degree polynomials. If we consider a 2 x 2 matrix, we may use M{z) = 
where ruj are assumed to be first-degree real (complex) polynomials in the real (complex) 
variable z and cTj are the Pauli’s matrices [23|. With this choice the condition M{zY = 
constant/ implies the following constraint 


[uAq\ du 


(17) 


Y(z) = y(0) 


/ -|- Mo / cosh [uAq] 
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3 

M{zf‘ = = {az^ -\-hz + c)I = constant/ ^ a = 0, 6 = 0. (19) 

i=l 

where a and 6 are functions of the coefficients of the polynomials rui. Since each 
polynomial has two coefficients we conclude that we have 2 x 3 — 2 = 4 degrees of 
freedom with respect to the choice of rui. Increasing the dimension of the matrix M it 
increases also the freedom’s degrees. This is due to the fact that conditions flTTll give 
two constraints, independently of the dimension of the matrix M. To better clarify 
this point, let us consider the matrix M{z) = where 7 " are the Dirac’s 

matrices [23]. We have 


3 

M{zY = + bz + c)I = constant/ ^ a = 0 , 6 = 0 .( 20 ) 

As before, conditions oa give two constraints so that the arbitrary coefficients 
are 2 x 4 — 2 = 6 . 

A more complex example. Let us consider the following conditions 


M{zf = 0, 


d 


= 0, M = —M{z), 


dz 


d^ 


M = constant, M =-—Miz). 

dz^ 

After a long but straightforward algebra we obtain 


( 21 ) 

( 22 ) 

(23) 


O^M{z) = TMhP, 

Oi”+‘M(s) = TMM'\ 

= 2” (m + mm) M". 


(24) 

(25) 

(26) 


00 


Mz) = > —- = —' ^0 “ +-cos 

3n! 3 

n =0 

00 3n+l 

Mz) S. = Af„ (3 X 


P/3. 


V3Al^\ 

22/3 


,(27) 


-1 


n=0 






u 


22/3 


7 / 3 .. 


+ vSe 22/3 


sm 


22/3 


/3(z)^5^03n+2^ 


u 


3n+2 


n=0 


e2^/^AV3._e-4^cos 


(3n + 2)! 

^A^u 


— ( il/n + il/nil/n ) ( 3 X 


U 


-1 


, (28) 


22/3 


y/3. 


- \/3e TW 


sm 




u 


22/3 


(29) 
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Y(z) = V(0) 



(30) 


Note that the integral of the f fnnctions can be easily performed. As before the solntion 
holds trne for an arbitrary dimension matrix M. 


4. Time dependent Hamiltonian 


The solntion of time dependent Hamiltonian is a fnndamental topic in qnantnm 
mechanics, qnantnm field theory, qnantnm theory of many particles etc. (see for example 
Ref. [1] as excellent text and review book). Using the previous results we may write 
a formal solution of the quantum problem under consideration. From Eq. ([5]) we infer 
that the operator O has the following expression 

O = A - (31) 

Using the formal solution fll2p . via Trotter formula, we may rewrite the exponential of 
the operator O as 


exp [uO] 


N 

lim I I exp 

V^oo 

i=i 




exp 




(32) 


where for simplicity we set h = 1 and it is understood that the operator exp acts 

on the left. Combining Eqs. flT^ and fl52]) we can write the ket corresponding to the 
solution of Eq. (|5]) as 


\A) =IV’o) + 


N 


SlU / |-*^2(io)lexp 


u 


exp 




|to=o du iV’o)- 


(33) 


t=i 0 

Note that the operator exp translates of a quantity ^ the argument of the 

Hamiltonian H{to). The wave function is 


fj{x, t) = {x\^jJt) = fj{x, 0) + 


N 


lim / II {x\[—iH (m)] exp 

N^oo I 


u 




u 

N 


du I'fo). 


(34) 


It is worthy to stress that in this formalism the time-dependence of the Hamiltonian, 
i.e. the dependence on to, is transformed in a parameter dependence. The integration 
variable that gives rise to the time dependence is u. 
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We achieve a first result. Taking the time derivative of Eq. (155]) . after little algebra, 
we have that for the state |'0t) holds 

N 

iV’t) = lim TT exp 

N^OO 

The above expression is nothing but that the solution to the Hamiltonian problem using 
a discrete approach 

We now assume that we know the orthonormal eigenstates of the Hamiltonian, 
i.e. Hlfo) |ntp) = Eniff) In^g) where, for what we said, the ’’time” fo has to be considered 
as a parameter of a time independent Hamiltonian. Using the closure relation [26] 


^ nf 

-i—H 3 


N 


N 


I'fo)- 


(35) 


K)(^iol= ("iioK) = ^nm, 

n 

we may rewrite Eq. 


as 


(36) 


tpix.f) = + 

N * 

lim > TT / du 


ni—riN i=l 


{nj+i\nj) [-zEnj^{u)]exp 


N 


U 


-t—J^En, [3 

j=0 


u 


^f{x,0), 


(37) 


where '0(a;, 0) = {x {tfo), and \nj) =\nj,jj^). Inserting a new set of variable |a;) with the 
properties 


(38) 


J dx\x){x\=I, {x\y) = 6{x — y), 

— OO 

we obtain the alternative expression 


'rp{x,f) = i>{x,0) + 

N ^ OO 

eh/*/ 


dxi 


ni—riN i=l 


{Xj) 'Ipn^ {Xj) [-lEnM] exp 


N 

iV ^ 


j=0 


N 


^f{x,0). (39) 


where f>nj{x) = ifuj {^^3jf) = (^1%)- Eqs. (1571) and (155]) represent the exact solution of 
the starting problem, Eq. ([5]). Making the approximation 




u \ u d 


U \ U 

nj+i,3j^), ^ = (40) 
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K+iK) 


. u 




nj+i,nj- 


( 41 ) 


The sums on the states | Uj) can be preformed using the delta Kronecker for two 
consecutive indices, and we have at hrst order approximation for N ^ oo 


t) ^ 0) + J2 u) [-^EM] 


exp 


-i J Eniz)dz\A 2 ) 

0 


where, by dehnition, c„ = {uq I-^q) and we used the continuous limit 


N “ 

j E„(z)dz. 


(43) 


i=o 0 

Via a simple integration by parts we obtain 


fj{x,t) Ri y^^Cnfjn{x,u) 


exp 


-i / En{z)dz 


+ 


t 

/ d 

'0n(a:,M)exp 


-i J En (z) dz 
0 


(44) 


The hrst term on the right side represents the adiabatic approximation (hrstly 
introduced by Born and Fock [27]) while the second is the correction to the adiabatic 
term. Let us now evaluate the next order approximation. Keeping only the terms 
containing the factor u/N we may write the correction as 


ISfj) = 



d 

dyEn{u)ck-^{n,z\k,z) x 


X exp 


Ek{p)dp+ / En{p)dp 


|n, u) 


(45) 


or in terms of wave function 


Sfi^x, t) 



d 

dyEn (u) Cktfkiy, z)—f)l{y, z) x 


X exp 


I j Ek (p) dp + J En (p) dp 
,0 2 : 


'ijjn{x,u). (46) 
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The integral on the space variable y can be expressed in terms of transition probability 
and we may write in more compact way 


d 


6'll){x,t)= / du / dzEn{u) Ck-^Pkn{z) X 

^ u j Jo 


n,k 


X exp 


-i\ Ek (p) dp+ En (p) dp 


iln{x,u). 


(47) 


OD 

where Pkn{u) = / 'if!f{x,u)'ipk{x,u)dx is the transition probability from k to n. Note 

— CX> 

that the diagonal terms k = n vanish identically. In order that the above expression 
represents a correction to Eq. (14^ it must hold 


II 5ii{x,f) ||<|| 'f{x,t) II (48) 

where 'ipfx^f) is given by Eq. fl42|) . The above inequality Exes the range of validity of 
our approach. As we will see in the next section, through this formalism, we can obtain 
the Wentzel-Kramers-Brillouin (WKB) approximation [TJ [23] for ordinary differential 
equations. 

5. Ordinary differential equations 

With some caution, the previous approach can be extended also to ordinary differential 
equations, namely when the matrix driven the system is not a hermitian operator. 

4r(() = M[t)Y(t). (49) 

In the previous section an approximated solution has been obtained using the closure 
relation and the orthonormal condition. Since the matrix of the system is not an 
hermitian matrix, we need to define two sets of vectors 


^ \n){n\= /, (m|n) = 6nm, (50) 

n 

where (rl| in general does not coincides with the hermitian conjugate of In). Repeating 
the same steps of the previous section, we end up in a similar approximated expression 


\Y,) « |K„) + ^ / duCnXn{u) exp 


An {z)dz 


\n{u)), 


(51) 


where Cn = (no | To) and An(n), |n(n)) are the eigenvalues and eigenvectors of M(n) 
respectively. For consistency with the previous section we also introduced the symbols 
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|Ft) =Y(t), |Yo) = ^(0). To better explain the procedure and the notation let consider 
the second order differential equation 

y = f{t)y. (52) 


Writing it in matrix formalism we have 


fL 

dt 


\Yt) 


0 1 
fit) 0 


lYt), I Yt) 



The eigenvectors and eigenvalues of the driven matrix are 


a 

b 


( 53 ) 


|l) = ]Vi 


, |2) = N2 


1 

h/W 

Ai = Vm, v = -VM 


-V7W 


( 54 ) 

( 55 ) 


with Ni a factor that has to be chosen in such a way to satisfy conditions (l50|) . It is 
straightforward to obtain 


iVl,2 = 
( 1 | = 


1 

TWm 

1 

7WW) 


Vm i], (2| 


Consequently we have 


1 

vWW) 


•/Tw -1 


( 56 ) 

( 57 ) 


yi{t) = y(t) ss o+ f 


/(o) 


X 


a\/ /(O) sinh 


\/ fiz)dz 


L^o 


+ b cosh 




uo 


du, (58) 


2 / 2 (t) = yit) ~ 5 + 




a\/ /(O)cosh 


^/Jiz)dz 


uo 


+ b sinh 


Vfi^)d< 


uo 


du. (59) 


Performing an integration by parts and redefining the constants f/ /(0)a = A 


and 


-1 


b = B we may write Eqs. fl58ll and fl59ll in the following form 


yit) 


A 


Vm 

1 f‘ /(«) 


cosh 


^/Jiz)dz 


uo 


B . , 

H- , sinh 


Vm 


4 


A cosh 


\fJi^dz 


uo 


+B sinh 


uo 

pU 


\/ fiz)di 

Vfi^)dz 


uo 


+ 

du, (60) 
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vit) 


A 


1 

4 


VM 

f{u 


sinh 




L~'o 


+ 


B 


VW) 


cosh 


Uo 


V 


[A sinh 


\/K^dz 


+B cosh 


\/ f{z)dz 

Vfi^)dz 


+ 


du. (61) 


Similar expression can be found in terms of trigonometric functions, according to the 
sign of f{t). Also, in the hrst term on the right side of Eq. (I60|) we can recognize the 
WKB approximation. We now consider a further example of differential equation and 
compare the exact solution with the approximated solution with the correction given by 
Eq. fl45|) . Making explicit formula fl45|l for the problem under consideration, we obtain 


\6Y{t))= lduX{u) / dz 




exp 


exp 


X{p) dp — I X (p) dp 


- I X{p)dp+ I X (p) dp 

Cl I 2) 


C2 I 1) + 
(62) 


where A(m) = /(«), f{z) = dzf{z), Cn = (no|E(0)) and the vectors In) are given by 

Eq. (1541) . Let us consider for example the following differential equation 

+ Vt + ly = 0, 2 /( 0 ) = 1, p(0) = 0. (63) 

We can compare the exact solution with the hrst order approximation given by Eq. (1581) 
and the approximation where we add the correction fl62l) . The result shows an excellent 
agreement when we add the correction (see Figs. [Hand[2]). 


y(t) 



Figure 1. Comparison be¬ 
tween the exact solution (con¬ 
tinuous line) and the approx¬ 
imated solution (dashed line) 
given by Eq. (1551) . 


y(t) 



Figure 2. Comparison be¬ 
tween the exact solution (con¬ 
tinuous line) and the approx¬ 
imated solution (dotted line) 
given by Eq. (1551) where it 
has been added the correction 
given by Eq. (1551) . 


More in general it can be shown that for a function asymptotically behaving as 
a power law, f(t) ~ ±t", with respect to the differential equation (1521) the following 
properties hold true (for sake of brevity we omit the analytical calculations): 
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i) In order that holds condition 0481) . for an asymptotic power law behaviour with 
f{t) > 0, it is sufficient that a > —2. This constraint on the power law coincides 
with the region of validity of WKB approximation.. 

ii) For f{t) < 0 approximated solution ([58]) gives the same asymptotic results of 
WKB approximation plus an extra constant for a > —2. For correction fl62l) 
to balance the extra constant it is sufficient that —2 < a < 1. In this range 
of the parameter a correction fl62l) improves the result with respect to WKB 
approximation. 

It is worthy to stress that condition (i) is sufficient but not necessary. A simple example 
is given by the discontinuous function f(t) = a9{l — t) where 6{z) is the step function. 
It is easy to check that formula (l58l) gives the exact solution. We can better appreciate 
the novelty of our approach studying problems in hnite interval containing the turning 
points, i.e. the points satisfying the equation fit) = 0. As we will see the accuracy 
of the presented approximation is very good and in the next section we will show an 
analytical formulation for the eigenvalue problem associated to a differential equation 
with boundary conditions. 

We will devoted the rest of the paper to this important problem. To better show 
this, we consider the following example where WKB approximation can not be used, at 
least directly. Let us study the following differential equation 


y(t) 



Figure 3. Comparison be¬ 
tween the exact solution (con¬ 
tinuous line) and the approx¬ 
imated solution (dashed line) 
given by Eq. (l5^ . 


y(t) 


5 - 

4 - ' 


3 



- 1.5 - 1.0 - 0.5 0.5 1.0 1.5 


Figure 4. Comparison be¬ 
tween the exact solution (con¬ 
tinuous line) and the approxi¬ 
mated solution (dashed line) 
given by the WKB approxi¬ 
mation. 


- {cos^ t)y = 0, t/(0) = 1, t/(0) = 0, te [-4,4]. (64) 

In Figs. [3] and m we see that the approximation given by Eq. fl58|) shows a small distortion 
with respect to the exact solution in the neighbourhood of the turning point t = n/2 
while WKB approximation, given by the formula 


vit) 




yj f{u)du 




(65) 
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except for a region near the origin, does not hold. As final remark of this section we 
note that the presented approach can be straightforwardly extended to higher order 
differential equations. 

6. Eigenvalue problem 

In this section we will face the problem of hnding the eigenvalues in the problem of a 
differential equation with boundary conditions. To show that, we will study in detail 
the following example of eigenvalue problem. Let us consider the differential equation 
dehned in the interval t E [— 1 , 1 ], that is to say 


df^ 


y + \{l- t^)y = 0, 2 /( 0 ) = 0, 2 /( 1 ) = 0, A > 0. 


The problem admit a general analytical solution given by 


( 66 ) 


y(t) = e ( aff, +bF J (l “ ^) . 5 . 


(67) 


where Hy{z) is the Hermite function and F{a,/3,z) is the Kummer confluent 
hypergeometric function. Imposing the condition 2 /( 0 ) = 0 we obtain 


y{t) = 


3 

4 


^ 3 

2 ’ 


'/Xt'^ 


( 68 ) 


while from the right boundary condition, i.e. 2 /( 1 ) = 0 , we obtain the numerical 
values of A that gives nontrivial solutions to Eq. fl 66 |) . In spite of the fact that WKB 
approximation is diverging at t = 1 still it gives hnite values for A that vanish the 
function at t = 1 


UWKBit) — 




sm 


\/ Xn 


(^tVl — f^ + arcsin t j 


, A„ = 16^2. (69) 


Note that the applicability of WKB approximation rests on the fact that the divergent 
factor is compensated by the function sinus. Alternatively, we may evaluate the 
eigenvalues using the approach presented in this paper. From Eq. fl58|) . via the 
condition 2 /( 0 ) = 0 , we have 


vit) 


^ a 

Jo 


cos 


^M\/l — + arcsin wj 


du. 


(70) 


Using the condition 2/(1) = 0 we have an equation for the eigenvalues. In Tabled] we 
compare the exact values of A with WKB and the present approximation. The percent 
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Xn (exact values) 

An given by Eq. flTOj) 

% error 

An given by WKB 

% error 

13.486 

13.767 

2 % 

16 

19% 

58.811 

59.174 

0 .6% 

64 

9% 

136.140 

136.557 

0.3% 

144 

6 % 

245.470 

245.930 

0 .2% 

256 

4% 

386.802 

387.296 

0 .1% 

400 

3% 


Table 1. A comparison between the exact values obtained from Eq. (l68l) and the 
approximated values given by Eq. (l5^ and WKB approximation. 


error is defined as | | where Xex is the exact value and Xappr is the approximated 

value. 

In the next example we will consider a case where, due to the divergence at the 
turning points, WKB approximation can not be used directly. Let us modify the problem 
with boundary conditions stated in Eq. fl66l) as follows 


+ A(1 - t^)y = 0, y{x) = y{-x), y{l) = y{0), A > 0. (71) 

Using the approximate expression for y{t) given by Eq. fl58|l . we may write the following 
equation for the eigenvalues 


/ (i-«b 


2,1/4 


sm 


(uVl — + arcsin u 


du = 0. 


(72) 


We obtain a sequence of of approximated eigenvalues in good agreement with the exact 
ones (see Table [2]). Note that the boundary condition is at the turning point, i.e. t = 1, 
and we can not use directly the WKB approximation since the approximated function 
is divergent at the turning point. 


An (exact values) 

An given by Eq. ([72]) 

Relative error 

43.185 

46.138 

7% 

77.736 

74.721 

4% 

208.573 

213.915 

2 .6% 

286.144 

281.010 

2 % 

500.880 

508.297 

1.5% 


Table 2. A comparison between the exact values An and the approximated values 
given by Eq. (17^ 


Finally we study a differential equation with mixed boundary conditions, i.e. 
Dirichlet-von Neumann conditions 


Ay + A(1 0, y(0) = 0, = j(l) = 0, A > 0. 


(73) 
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Evaluating the solution at t = 0 we obtain again expression 0681) . Vanishing the 
derivative of y{t) at t = 1 we obtain the eigenvalues of the problem. Since the condition 
on the derivative is at the turning point, i.e. t = 1, we can not use directly the WKB 
approximation. On the other hand we may use the approximate expression for the 
derivative of y{t) given by Eq. fl5^ that yields 


y{t) ^ b 



2 



+ arcsin 



• (74) 


The values of A so obtained are compared with the exact ones in Table [3l In general 
we can use Eqs. fl5^ and fl5^ as formulas to hnd, analytically, the eigenvalues of a 
differential equation with boundary conditions defined in a finite interval. 


An (exact values) 

An given by Eq. fl74|) 

Relative error 

5.122 

4.721 

8 % 

39.661 

39.836 

0.4% 

106.249 

106.063 

0 .2% 

204.856 

204.952 

0.05% 

335.473 

335.352 

0.03% 


Table 3. A comparison between the exact values obtained vanishing the derivative of 
Eq. (1551) at t = 1 and the approximated values given by Eq. dZlD 


7. Concluding Remarks 

In this paper it has been showed a formal analytical solution for a linear differential 
equation of arbitrary order n and with variable coefficients. The solution has been found 
exploiting the reduction of the problem to a first order differential equation through 
the matrix formalism. This approach provided a solution also for differential equation 
involving operators such as the time-dependent Schrodinger equation. 

It has been shown that the most important approximations available in literature 
can be deduced from the proposed solution. With respect to ordinary differential 
equation, a noteworthy fact is that the proposed approximation is smooth around the 
turning points and this allows to study solutions of differential equation at points where 
the WKB approximation, in general, can not be used. Examples of analytical eigenvalue 
evaluation have been showed. Even though we studied examples of second order 
differential equations, the approach can be easily extended to higher order differential 
equations. Finally this work showed that there are several possible research lines 
that have to be fully explored. In particular the approach based on Baker-Campbell- 
Hausdorff formula is left for a future work. 
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